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Vacuum type of SU(2) gluodynamics in maximally Abelian and Landau gauges 
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The vacuum type of SU (2) gluodynamics is studied using Monte-Carlo simulations in maximally 
Abelian (MA) gauge and in Landau (LA) gauge, where the dual Meissner effect is observed to work. 
The dual Meissner effect is characterized by the coherence and the penetration lengths. Correla- 
tions between Wilson loops and electric fields are evaluated in order to measure the penetration 
length in both gauges. The coherence length is shown to be fixed in the MA gauge from measure- 
ments of the monopole density around the static quark-antiquark pair, ft is also shown numerically 
that a dimension 2 gluon operator A~ {s) and the monopole density has a strong correlation 
as suggested theoretically. Such a correlation is observed also between the monopole density and 
A'^(s) = A~{s) -\- A^^s) condensate if the remaining (7(1) gauge degree of freedom is fixed to 
(7(1) Landau gauge (UlLA). The coherence length is determined numerically also from correlations 
between Wilson loops and A^ A~ {s) and A^{s) in MA -I- UlLA gauge. Assuming that the same 
physics works in the LA gauge, we determine the coherence length from correlations between Wilson 
loops and ^{s). Penetration lengths and coherence lengths in the two gauges are almost the same. 
The vacuum type of the confinement phase in both gauges is near to the border between the type 
1 and the type 2 dual superconductors. 

PACS numbers: 12.38.AW,14.80.Hv 



INTRODUCTION 



It is conjectured that the dual Meissner effect is the 
color confinement mechanism P, The conjecture 

seems to be realized if we perform Abelian projection Q 
in the maximally Abelian (MA) gauge |3, • Abelian 
component of the gluon field and Abelian monopoles 
are found to be dominant 0, Q- Abelian electric field 
is squeezed by solenoidal monopole currents 0, ITol| . 
Monopole condensation is confirmed by the energy- 
entropy balance of the monopole trajectories [ll| and by 
evaluation of the monopole creation operator AH 
these facts support the conjecture that the color con- 
finement is due to the dual Meissner effect caused by 
the monopole condensation. Numerical calculations show 
that the vacuum of quenched SU{2) QCD {SU{2) gluo- 
dynamics) is near the border between th e ty pe 1 and 
the type 2 dual superconductor [l^, HI; Hll i although 
there are some claims that it is a superconductor of type 
1, Ref. H il6.] . Since the definition of a dual Higgs field is 
unknown the coherence length was calculated using clas- 
sical Ginzburg-Landau equations, while the penetration 
length can be calculated directly measuring the correla- 
tions between Wilson loops and Abelian or non- Abelian 
electric fields. 

In this note, we show that the coherence length 
could be derived directly from the measurement of the 
monopole density around a chromomagnetic fiux spanned 
between a static quark-antiquark pair. We use the dual 
Ginzburg-Landau effective theory of infrared SU (2) gluo- 
dynamics |l7l , evaluate the monopole density around 



the flux theoretically and compare it with the value ob- 
tained numerically. 

We consider also the dimension 2 gluon operator 

A+A-{s)^Y.^{Al{s)f + {Al{s)f] 

in the MA gauge. The MA gauge is a gauge which 
minimizes a functional '^gA'^A~{s). It is well known 
that the off-diagonal gluon fields A^^{s) with z = 1, 2 are 
small everywhere except at the sites where monopoles 
exist. Hence a strong correlation between A'^A~{s) and 
monopole currents |fcp(s)| is expected. The off-diagonal 
gluons have no essential effects on the confinement of 
fundamental charge, whereas they can explain the screen- 
ing of adjoint charge If we perform the additional 
UlLA gauge fixing after the MA gauge is fixed, the op- 
erator can have a physical meaning. It is 
expected f^rom the previous results suggesting monopole 
dominance^ that monopole contribution could explain 
all non-perturbative effects in the quantity 



Hence we expect that the coherence length can also be 
derived from correlations between the Wilson loops and 
the dimension 2 operator A+A~(s) or between the Wil- 
son loops and the dimension 2 operator 

A^{s)^A+A-{s)+A^A^{s). 

We find that the coherence lengths determined from the 
monopole density and the dimension 2 operators are con- 
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sistent with each other. We also observe that the pen- 
etration length and the coherence length are almost the 
same and we conclude that the vacuum is near the border 
between the type 1 and the type 2 dual superconductors 
in the MA gauge. 

Since the MA gauge - in which the confinement mech- 
anism is definitely found to be realized - is just one gauge 
among infinite possible gauges and, on the other hand, 
the physics should be gauge-independent it is important 
to know the confinement mechanism as well as the type 
of the vacuum in another gauge and in the case in which 
the gauge fixing is not performed at all. 

This problem has been discussed recently in Ref. 0| 
where the Landau (LA) gauge is considered and for 
Abelian components the dual Meissner effect is observed. 
A magnetic displacement current plays the role of the 
solenoidal super current which squeezes the Abelian elec- 
tric fields, although the density of DeGrand-Toussaint 
monopoles [IJ is very small in the LA gauge. The ob- 
servation of the dual Meissner effect in the LA gauge 
suggests that there exists a gauge-independent definition 
of the monopole and, consequently, of the monopole con- 
densation. There are some attempts to find a gauge- 
independent definition of magnetic monopoles [23, 123, 
|2^ . Based on the analogy of the SU(2) gluodynamics in 
the LA gauge with a nematic crystal in Ref. 25] an ex- 
istence of various topological defects was suggested. But 
the definite answer about degrees of freedom which are 
relevant to the confinement in the LA gauge is not yet 
obtained. 

It has been shown that non-perturbative part of the 
condensate (A^(s)) is explained completely in terms of 
monopoles in compact QED in Landau gauge |2^ , where 
the monopole condensation is known to be responsible 
for the confinement of charge . The non-perturbative 
part of {A'j^{s)) corresponds just to the vacuum expec- 
tation value of a dual Higgs (monopole) field. For glu- 
odynamics the relevance of the M^ (s)) condensate for 
confinement is discussed in Refs. "2^, '29'|. Using the sum 
rule technique the mass gap generation due to d = 2 
gluon condensate is discussed in Ref. [30l| . 

In this note we fix the type of the vacuum also in the 
LA gauge. First we measure the penetration length from 
the electric field flux as done in the MA gauge. Then 
we fix the coherence length from correlations between 
Wilson loops and the dimension 2 gluon operator A'^{s), 
assuming that a relation between the dimension 2 opera- 
tor and an unknown gauge-independent monopole exists 
in the LA gauge similarly to the MA gauge. We find 
both the penetration length and the coherence length in 
the LA gauge are consistent with those in the MA gauge. 
The type of the vacuum is found to be gauge- independent 
as it should be. Note that the LA gauge corresponds to 
a gauge in which the functional ^^(s) is minimized. 
Thus the operator ^^(s) could have a physical mean- 
ing in LA gauge p9l Isij if the Gribov-copy problem is 
solved. 



II. CONSIDERATION IN THE DUAL 
GINZBURG-LANDAU THEORY 

A. General dual Ginzburg Landau picture 

The monopole density around the QCD string is 
described very well by the dual superconductor pic- 
ture 0, |32- The dual superconductor (or, the dual 
Ginzburg-Landau (DGL)) Lagrangian corresponding to 
SU(2) gluodynamics has the following form: 

+<m'-v')\ (1) 

where i?^ is the dual gauge field with the mass ms — grj, 
and $ is the monopole field with magnetic charge g and 
with the mass m$ = a/Sk?]. In the confinement phase 
of SU(2) gluodynamics the monopoles are condensed, 
I < $ > I = ry. The coupling k defines the quartic inter- 
action of the monopole field $. Below we discuss some 
general well-known properties of the Abrikosov-Nielsen- 
Olesen [s^ vortex in this Abelian model. 

There are two mass-scales in the discussed Abelian 
Higgs model: the coherence length ^ and the penetration 
length A, which are inversely proportional to the masses 
of the monopole and the dual gauge boson, respectively: 

, A = — . (2) 

The border between the type 1 and type 2 superconduc- 
tors is defined as a region in the phase diagram space 
where the both length coincide, £, = X- 

We are interested in the behavior of the monopoles 
around the QCD string. The classical equations of mo- 
tion of the DGL model |^ contain a solution correspond- 
ing to the QCD string with a quark and an anti-quark at 
its ends. The infinitely separated quark and anti-quark 
correspond to an axially symmetric solution of the string. 
We consider the static solution which is parallel to the 
third direction of the reference system, 

Hp) - vfip)e'\ 

B, = ^^h{p), 53 = 0, B4 = 0, (3) 
9 P 

where f{p) and g{p) are dimensionless functions of the 
transverse distance p = y/x\+x^ to the string, and e.^ is 
the standard fully anti-symmetric tensor, = —e-ji and 
£12 ~ 1. The azimuthal angle is[4^ = arg(a;i -f ix2). 
Both functions / and h of Eq. ^ tend to zero as p — > 
and they approach the unity as p oo. 
The DGL classical equations of motion are: 

DliB)^ - A^m^ - V^)^ = , (4) 

[dlS^, ~ d^d,)B, = gk^i<P,B) , (5) 

where fc^ is the monopole current given by the following 
expression: 

= 3m[$*i?^(B)$] = m\d^ arg$ + gB^) , (6) 
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where D^{B) = + ig is the covariant derivative. 

In terms of the functions / and h used in the ansatz ||2J), 
the current © is given by: 



fe=0, fc4-0.(7) 



P P 



To derive this equation one should use the relation 
dO/dxi = —eijXj/p^. 

In terms of the ansatz Q the classical equations of 
motion are: 

np) + ^-#[i-Mp)F 



+ ^[1-/'(P)]/(P) = 0, (8) 

h'\p)-^+ml[l-h{p)]f{p)^Q. (9) 
P 

Expanding the profile functions at large p, f{p) = 1 — 
Sf{p) and h{p) = 1 — Sh(p), and keeping only linear 
terms in Eq. ^ and Eq. (|SJ|, we get the linearized 
classical equations of motion: 



Sfip) 



^-m|.<5/(p) = 0, (10) 



Sh"{p) - - ml ■ 6h{p) = , (11) 

P 

which have the solutions: 

^f{p) = CfKo{mBp) , Sh{p) = Chm^pKi{mq,p) .(12) 

Here Kn are the modified Bessel functions with the fol- 
lowing asymptotic (x — > oo) expansion: 



l + 0( 



(13) 



For the string solution with a lowest non-trivial flux the 
arbitrary coefficient C/ is always equal to unity, C/ = 1, 
while the coefficient Ch is equal to unity in the Bogo- 
mol'ny limit (i.e., exactly on the border between the type 
1 and type 2 superconductivity) , Ref. ,3^ ,35} . Since the 
numerical results suggest strongly that the SU(2) gauge 
theory is close to the border, we set Ch = 1 in our qual- 
itative discussion below. Thus, the functions h and / at 
large values of p behave as follows: 



fip) = 1 - loim^p) + ■■■ 
h{p) = 1 - niBphimBp) 



(14) 
(15) 



The QCD string is well described by the solutions 
of the classical equations of motion of Lagrangian 
The qualitative behavior of the monopole condensate, the 
electric field and the angle component ke of the monopole 
current around the QCD string are shown in Fig. ^a), 
(b) and (c), respectively. 

Summarizing, the value of the monopole current near 
the QCD string (obtained from the classical equations 



of motion) is zero in the center of the string and it is 
also zero far from the string. The current has a maxi- 
mum at a certain distance (which is numerically found 
to be about 0.2 fm in the DLG corresponding to SU(2) 
gluodynamics [l^ ^2J). The only non-zero component 
of the classical monopole current is the angle compo- 
nent kg, while other components (radial, temporal and 
z-component) are zero, fcr- = 0, A;4 = and ^3 = 0. 



B. Monopole density around QCD string 

In the numerical calculations the distributions of the 
monopole current around the QCD string E is measured 
with the help of the correlation function 



^string 



{k^iO)Wc) 
{Wc) 



as = c, (16) 



where S denotes the string world-sheet corresponding 
to the minimal surface spanned on the Wilson contour 
C. The expectation value (|16(l is non-zero contrary to 
Eg . (|47|l due to the broken Lorenz invariance because of 
the presence of the string. 

The monopole density is non-zero in the absence of 
the string. We call this value of the density as "vacuum 
monopole density", There are two contributions 

to this monopole density coming from (i) the long (in- 
frared) monopole loop which forms the monopole conden- 
sate |2a,E3 and from (ii) the small monopole loops which 
represent the perturbative (ultraviolet) fluctuations. 

Very naively, the presence of the string should make 
the monopole density bigger: the vacuum contribution 
gets an additional contribution coming from the classical 
(solenoidal) current k"^^^^^ = k^*^"'^^. The naive picture is 
plotted in Fig. E| 

Thus, naively, the density of the monopoles should 
increase at some distance from the string. Moreover, 
naively one expects that at large transverse distance 
p from the string the monopole density, [according to 
Eqs. H7I13I14I15|| ] is controlled by the penetration length 
since ^ |fcvac| _|_ const exp{-TOBp}. 

However, the described qualitative picture definitely 
contradicts the numerical results obtained in 

Ref. [33. In order to investigate the behavior of the 
monopole density near the QCD string we study analyt- 
ically the London limit in the next subsection. 



C. Monopole density in the vicinity of QCD string 
in the London limit 



The London limit is characterized by the infinitely 
deep potential, m$ — > 00. The Lagrangian of the DGL 
model in this case is 



^dgT = \ (% S.])' + y (5^^ ^gB^r . (17) 
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Classical 
solution of DAHM 




classical 
solution of DAHM 



1/m 



(a) 




? 7/^ 



FIG. 1; The qualitative behavior of (a) the monopole condensate "I>, (b) the electric field Ez and (c) the magnetic current 
component ke around the center of the string as functions of the transverse distance p. The string is given by a classical solution 
of the DGL model. 







Naive picture 




FIG. 2: The naive qualitative behavior of the density of 
monopoles around the QCD string. 



The QCD string S manifests itself as a singularity in the 
phase of the Higgs field: 



(18) 



where the string S is parameterized by the vector 
Xfi,{Ti,T2) which depends on the parameters ri^2- The 
antisymmetric string current is given by the follow- 
ing expression: 



The partition function of the model p7ll can be rewritten 
as a sum over the closed strings 4S] [39lj : 

Z = r Dip r HB exp{- / d4/:L™i™(i?,^)} 

= / DScxpj-S'.trlS)}, (20) 
where Sstr is the string action: 

^str(S) = 

27tW fd^x I d4yE^,(a;)i?„,^(x-y)E^,(y)(21) 



and D]\[ is the propagator of the massive scalar par- 
ticle, (-0)2 + M^)Dm{x) = (5W(a;). The strings are 



closed: d^'S, 



0. The derivation of the right hand 



side in Eq. H2U|I is easy to make by fixing the uni- 
tary gauge, if = and, consequently, making the shift 
Bf^ — > B|^^ — {l/g)dft(p. Then Eq. lfTH|) implies that un- 
der the shift dif,B^] d^^B^^ - (27r/g)*I]^^. Finally, 
having integrated over the Gaussian field B^ we get the 
right hand side in Eq. 

The sources of the electric flux (i.e., the quarks) run- 
ning along the trajectory C are introduced with the help 
of the Wilson loops written in terms of the original gauge 
fields A^. The quantum average of the Wilson loop Wc 
can be rewritten as a sum over the strings similarly to 
Eq. 



{Wc) = 

4 / DS exp{-5,tr(S) - ^currc„t(/)} , (22) 

where the strings are spanned on the current jc '■ S^S^i/ = 
j^. The action for the currents is given by the short- 
ranged exchange of the dual gauge boson: 

>5'currcnt(jc) = 

i fd^x [ d'yfJx)Dr., ix - y)£{y) , (23) 
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where e = 27r/g is the electric charge. 

Below we evaluate the density of the monopole current 
in the vicinity of the fixed QCD string. To this end we 
assume that the leading contribution of the QCD string 
is naturally given by the minimal surface configuration. 
Moreover, to avoid boundary (quark) effects, we place 
the static quarks at the (spatial) infinities of the axis x^. 
Consequently, the quark term (|23|l does not play any role 
in the forthcoming discussion. 

Thus, we consider the infinite static string which is 
placed along the third direction. The corresponding 
string current - calculated from Eq. 119|l - is given by 
the equation: 

Spy = (^(5p,3'5/j,4 - S^^3S^^i^d{xi)S{x2) ■ (24) 

The monopole current JHl in the London limit is 

k^^if[d^^ + gB^). (25) 
Let us consider the following generating functional: 

/oo 
-OO 



X exp 



L 



London 
DGL 



(S,^s) 



},(26) 



where the singular phase (^e corresponds to the string 
position S fixed via Eq. I|18|l . The monopole current 
in the presence of the string is given by the variational 
derivative: 



1 



Analogously, the (squared) monopole density is 



{x). 



c=o 



(27) 



(28) 



Proceeding similarly to the derivation along 
Eas. l2UI22|) wc get the following expression for the 
generating functional: 



z[j:,c] 

expj-^ < 



C^ix)D^f{x-y)CAy) 



2mri^C^,{x)D^f [x - y)a„*S„,(y) 



-5. 



tr(S)}(29) 



where D^j^{x) is the propagator of the massive vector 
boson Bfj,, and the string action is given in Eq. (|21|l . 

An evaluation of the vacuum expectation value of the 
monopole density (|27|) gives: 



I string 



-2^Tf I d^yD":j>{x-v)d^*i:^M.m 



In particular, in the case of the static string (|24() we get 
the classical London solution: 



strnig 



.string 



, fcf "'^s = , 



(31) 



where 



(32) 



is the propagator for a scalar massive particle in two 
space-time dimensions. Using Eas. (l31l32|l we get the 
explicit form of the only non-zero component of the 
solenoidal current: 



^string ^ rj'^rnsKiimBp) 



(33) 



The monopoles form a solenoidal current which circulates 
around the string in transverse directions. 
The squared monopole density is: 



where 



(fcu)s = (Ml + (fcu)o, 



(34) 



2 4 A 2 

167r2 



0(log(A/mB)) (35) 



is the quantum vacuum correction. We have regularized 
the divergent expression of the vacuum correction by the 
momentum cutoff A. The correction is quadratically di- 
vergent. 

The total (squared) density of the monopole current is 
given by 



rj^MlKfiniBp) 



167r2 



0(log(A/mB)) , (36) 



where the solenoidal current k^'^™^ in the London limit is 
given by Eq. (|31|l . This expression is the exact in the Lon- 
don limit (up to logarithmically divergent but distance- 
independent corrections). 

One may easily see from Eq. H3()|l that the naive ex- 
pectation of the density behavior - shown in Fig. [21- is, 
in fact, correct in the London limit. Then the total den- 
sity (in which the coherence length is zero) must have a 
maximum at the distance of the order of the penetration 
length, l/niB- However, the naive picture depicted in 
Fig. [5] is not valid the case of the finite coherence length 
considered below. 



D. Monopole density in the vicinity of QCD string 
for finite coherence length 

Here we show, that in the real case os a finite coherence 
length, the naive picture, described in the previous Sub- 
section, is no more correct. Indeed, in this case the value 
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of the monopole condensate is varying in the vicinity of 
the string, and the (quahtative, at least) generahzation 
of Eq. should read as follows: 



(fc; 



)2 + 

^2 . 9'mp)\'^' 



(37) 



where we have taken into account the behavior of the 
monopole condensate by the simple replacement |49| r] — > 
|$(p)| in 



167r2 



(38) 



Note that the quantum correction to the squared 
monopole current in the vicinity of the string (with p ^ ^) 
is not equal to the vacuum expectation value measured 
far outside the string (p » ^)! 

The quantum correction is much stronger than the 
classical one, therefore the leading behavior of the total 
density is controlled by the quantum corrections. The 
behavior of the monopole density in the vicinity of the 
string is shown in Fig. |21 by the solid line. The various 



Real picture 




FIG. 3: The real qualitative behavior of the density of 
monopoles around the QCD string. 

contributions to the total density are also shown in this 
Figure (the dashed lines). The theoretical expectation 
- shown in Figj_pi - is in agreement with the numerical 
result of Ref. "3^ 

Thus, we expect that the quantum corrections play 
an essential role in the case of a finite coherence length. 
Moreover, the leading behavior of the monopole density 
at large distances is controlled by the coherence length ^ 
(and not by the penetration length A). This fact can be 
seen from Eqs. (|13I14I15I) in the limit p > 



{\k\hip) 




(39) 



As it is discussed in Section IIIII the monopole den- 
sity should locally be correlated with the condensate 
(one can naturally expect a correlation with a 
short length scale of the order of 0.05fm - this topic will 
be discussed in the next Section). Therefore, the corre- 
lation of the monopole density with the QCD string H39(l 
indicates that the A~ condensate is also correlated 
with the QCD string. The correlation lengths for the 
"monopole density" -" string" correlations and for the 
" A^A~ condensate" -" string" correlations should be the 
same and equal to the coherence length of the dual su- 
perconductor: 



{A+{p)A-{p)h 



~ 1 — const • e 



(A+{p)Aj.ip))o 
with p 3> This is the main result of this Section. 



(40) 



III. A+A- CONDENSATE AND ABELIAN 
MONOPOLE IN MA GAUGE 

The MA gauge is defined by the maximization of the 
functional 



R[U]=yRi[U], 



I 



Uia^UJaa 



with respect to the gauge transformations. 



max R\U' 



Ir^ — fl^TI ft 1 - 



(41) 



(42) 



Using the standard parameterization of the SU{2) link 
matrices. 



cos (pi e 



sm (pi e 



(43) 



we obtain Ri[U] = cos 20;. The maximization makes (p 
as close to zero as possible. 

The off-diagonal fields, Uj^ = ± sin (pi e^'-^' correspond 
to continuum fields ±A^{x) and continuum quantity 
A^{x)A~{x) corresponds to the lattice quantity sin^ = 
{1 — Ri[U])/2. Thus we are able to make an identification 
(no summation over p is assumed): 



A+{x)A-{x) = -{l-RxAU])^ 



(44) 



where the equality is exact in the naive continuum limit. 

The first measurements of the local correlation between 
monopoles and the quantity Ri were done in Ref. [43l| . 
where the quantity Rc was used: 



Rc[U] = Y,Mu]^ 

ledc 



(45) 



The summation is going over all links belonging to the 
cube c. The distribution of the quantity Rc at the cubes 
occupied by monopoles and the cubes, not occupied by 
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monopoles, is observed. It is shown that at the monopole 
position the quantity Rc is generally smaller compared 
to the same quantity in the empty space. Therefore, 
the monopoles suppress the quantity i?c, and, accord- 
ing to Eq. H44() . the A~ condensate is enhanced on 
monopoles. 

One can suggest, that the correlation of the A^A^ 
condensate with the monopole is short ranged. Indeed, 
the correlation of the monopole with the SU(2) action 
in the MA gauge is short ranged, with the character- 
istic correlation length Q CAction ~ 0.05fm. Since 
the SU(2) action involves the off-diagonal components, 
it seems natural to suggest that the correlation length 
Ccond of the off-diagonal components of the gluon field 
(or, of the A'^ A~ condensate) with the monopoles is 
not much higher than the CAction- Thus, one can expect 

that Ccond ~ CAction ~ 0.05fm. 

Thus the A+^^-monopole density correlation func- 
tion. 



C(r) 



{\KmK^r)A-{r)) 
{\Km){At{Q)A-,{Q)) 



(46) 



at large r is an exponentially decaying function with char- 
acteristic length scale C- 

In general there are two types of correlations: along the 
monopole current fc^ and perpendicular to the monopole 
current. In Eq. H4t)|) we assume that the distance r is 
perpendicular to the direction of the monopole current, 
R^k^j, = {i.e., the correlations are studied in the trans- 
verse to the monopole current directions). Obviously, due 
to the scalar nature of the A~^A~ operator the correlation 
of this quantity with the monopole current is zero: 



B. The MA gauge case 

Non-Abelian electric fields are defined from 1x1 pla- 
quette Uf„y{s) = U°^{s) + iU^^{s)a°- as done in Ref. g||: 



fc)+C/4\.(s)] 



(48) 



The static quarks are represented by the Wilson loop 
W{R,T). The measurements of the electric field are 
mainly done on the perpendicular plane at the midpoint 
between the quark pair. A typical example is shown in 
Fig. ^ Note that electric fields perpendicular to the QQ 
axis are found to be negligible. 
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r[fm] 

FIG. 4: The non-Abelian E electric field profile in the MA -|- 
UlLA gauge obtained with the use of the J? x T = 6 x 5 Wilson 
loop. The solid line denotes the best exponential fit by the 
function 1491 1 with the best fit parameters given in Table |T| 



IV. NUMERICAL RESULTS 
A. Method 

We use an improved gluonic action found by 
Iwasaki ^3 which was already implemented in Ref. |20j |: 

5-^/3 |Co^ Tr(plaquette) -I- Ci^ Tr (rectangular) | . 

The mixing parameters are fixed as Co + 8Ci ~ 1 
and Ci = —0.331. We adopt the coupling constant 
/? = 1.2 which corresponds to the lattice distance a{(3 = 
1.2) = 0.0792(2)fm. The lattice size is 32* and we 
use around 5000 thermalized configurations for measure- 
ments. To get a good signal-to-noise ratio, the APE 
smearing technique is used when evaluating Wilson 
loops W{R,T) = W° + iWa". The thermahzed vac- 
uum configurations are gauge-transformed in the MA(-t- 
UlLA) gauge and in the LA gauge. In the LA gauge 
the functional Tr[Ufj_{s) + Uj^{s)] is maximized with 
respect to all gauge transformations. 



The correlation length is determined by an exponen- 
tial fit of the electric field 021 for large r regions, where 
r is a distance perpendicular to the QQ axis. Below we 
observe that the electric field as well as other field distri- 
butions around the string can be fitted well by a simple 
exponential function 

/(r) = ci exp(-r/C) + ca , (49) 

where C and ci_2 are the fit parameters. The correspond- 
ing best fit parameters are presented in Table ^ The 
best fitting curve for the distribution of the electric field 
is plotted in Fig. Bias solid line. From this fit we fix the 
penetration length 5u\ X. 

We show the results for the penetration length in the 
MA + UlLA gauge in Fig.|31for various sizes of Wilson 
loop in space directions R. Here we see the penetration 
lengths for both Abelian Ea and non-Abelian E electric 
fields are compatible with each other. This is expected, 
since in MA gauge off-diagonal gluon components are 
forced to become as small as possible. 

Next we study the correlation between the monopole 
density |fc^(s)| and the operator A~^A~{s) given by 
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Quantity 


gauge 


Fig. 


C [H 


Cl 


C2 


x'/d.o.f. 




MA + UlLA 


lil 


0.140(3) 


0.044(2) 





3.06 


i\k^\{A+ A-)u) 


MA 


IE] 


0.0606(9) 


1.08(3) 


-0.0103(2) 


0.003 


{Wk^) 


MA 


m 


0.10(1) 


0.016(7) 


0.014300(8) 


1.35 


{W(A+A-)u) 


MA 


m 


0.094(8) 


0.04(1) 


0.40967(2) 


0.006 


{W{A+A-)e) 


MA + UlLA 


uni 


0.106(9) 


0.12(4) 


0.49068(4) 


0.01 




LA 


m 


0.139(1) 


0.038(1) 





1.88 


{WAl) 


LA 




0.118(4) 


0.09(1) 


0.74023(4) 


0.003 



TABLE L The best fit parameters corresponding to the fits of various quantities by the function 1)49^ . We indicate the gauge 
where the quantity is calculated and the figure number where the quantity is plotted. We set C2 = when the best fit value of 
C2 is consistent with zero. 
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FIG. 5: The penetration lengths for non-Abelian E and 
Abehan Ea{0 and 9^^) electric field profiles in the MA + 
UlLA gauge and their i?-dependence. Here 9 is defined by 
[/p(s) = exp(i6i;j(s)cr°) and 6''*^ is given in Eq. 
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FIG. 6: The correlation between the monopole density |fe^| 
and the operator A'^ A~ in the MA gauge. The solid line 
denotes the best exponential fit by the function l)49|l with the 
best fit parameters given in Table H] 



Eq. iBH) . The correlation data is plotted in Fig. It 
is completely consistent with the theoretical expectation 
discussed in Section UTTl In particular, the scale of cor- 
relations between |fcu(s)| and A'^A~{s) is about 0.06 fm 
according to Table ^ This value is pretty close to the 
scale CAction ~ O.OSfm of the monopole-action correla- 
tions. 

The histograms of the quantity i?c, Ea. (|^ . are plot- 
ted in Fig. [3 We discriminate between the histograms 
obtained at the cubes unoccupied by monopoles and ob- 
tained at the cubes occupied by long infrared monopoles. 
From this figure, we can clearly see the enhancement of 
the A~^A~ condensate on the Abelian monopoles. 

Let us next derive the coherence length in the MA 
gauge. The correlations between the Wilson loop and 
the monopole density, and between the Wilson loop and 
the quantity A^A~{s) are plotted in Fig.|Sland in Fig. El 
respectively. 

In this gauge, the quantity A~^A~ (s) may have a physi- 
cal meaning as we have discussed above. If the remaining 
f/(l) symmetry is gauge fixed by the U{1) Landau gauge, 
the dimension 2 gluon operator A'^(s) also acquires a 



physical meaning. 

In order to study different definitions of the quantity 
A'^A~{s), we plot the profile of the A+A^(s) condensate 
in Fig. 1101 using another definition 

A+A, =J2{{[elis)r + [elis)r}) , (50) 

which uses the angles 0fx{s) given by the relation U^(s) = 
exp(«0°(s)cr'*). In Fig. we show the coherence length 
determined by the use of the quantities A+A^, A'^Ag 
and k^. From Fig. ^2 we conclude that within the error 
bars these coherence lengths coincide. 

It is interesting to determine a non-perturbative con- 
tent of the gluon operator. To this end we measure only 
the monopole contributions to the dimension 2 gluon op- 
erator A'^{s) after the MA gauge, and, subsequently, the 
additional U{1) Landau gauge are fixed. This way of 
defining of the non-perturbative quantities is justified be- 
cause it is known that in the MA gauge the monopole con- 
tributions are responsible for essential non-perturbative 
physics. The monopole contribution to the coherence 
lengths is plotted in Fig. E| 
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side the monopole locations. The vertical line indicates the 
vacuum average (-Rc). 
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FIG. 8: The correlation between the 7? x T = 6 x 5 Wilson 
loop and the square monopole density in the MA gauge. The 
solid line denotes the best exponential fit by the function (I49II 
with the best fit parameters given in Table 
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the 7? X T = 6 X 7 Wilson loop and the A^AZ condensate lHH . 
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FIG. 10: The same as in Fig.|Hlbut for the correlation between 
the _R X T = 6 X 7 Wilson loop and the A'^Ag condensate lICT 
in the MA + UlLA gauge. 



quantum-mechanically. 



We plot the dependence of the coherence lengths of the 
spatial extension R of the Wilson loop in Fig. One 
can clearly see that these values are almost independent 
on R. Moreover, it is very interesting that the values of 
the coherence lengths are almost the same as those of 
the penetration lengths. Hence, if this relation holds for 
larger R - which is very plausible - then the type of the 
vacuum must be near the border to the type 1 and type 2. 
This is consistent with the results of Ref. [Hill 11113 
obtained with the help of the classical equations of mo- 
tions of the dual Ginzburg-Landau theory . However 
it should be emphasized that in this paper the determi- 
nation of the vacuum type was done for the first time 



C. The LA gauge case 

The Abelian electric fields E"^^ are defined in the LA 
gauge similarly to the MA gauge. We use the Abelian 
plaquettes d'^^(s) defined with the help of the link vari- 
ables 9'^{s): 

o%{s) ^ eiis) + ei{s + A) - ei{s + 6i{s) (51) 

where 6*^(5) is given by [/^(s) — exp(i6' °(s )cr°). 
The dimension 2 gluonic operator is |51|: 

A-{s)^Y.Y.{ei{s)f . (52) 

11— 1 a— 1 
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FIG. 11: The coherence lengths ^' for A+A:^ and A+Ag and 
k"^ in the MA + UlLA gauge. 
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FIG. 12: The coherence lengths of dimension 2 gluon opera- 
tor in the MA + UlLA gauge for various R. 



First let us discuss a determination of the electric fields 
around a pair of a static quark and an antiquark in the 
LA gauge. Since the confining behavior of the chromo- 
electric string is seen for large enough quark-antiquark 
distances R, we have performed the measurements for 
various R and T. A typical example is shown in Fig. 1141 

The R dependence of the penetration lengths is shown 
in Fig. 1151 The maximum quark distance in Fig. 1151 is 
R = 0.71fm which may not be large enough to see the 
confining string behavior. On the other hand, we see 
a small but clear discrepancy between the penetration 
lengths of the Abelian Ea and the non- Abelian E electric 
fields. The authors think it is caused by too small dis- 
tance between the quark and anti-quark so that the dif- 
ferent effects from the Columbic interaction may still play 
a role. Anyway for the interquark distances R < 0.71 fm, 
the Abelian dominance is not observed in the LA gauge. 
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FIG. 13: The coherence lengths of dimension 2 gluon op- 
erator and the penetration length in MA -j- UlLA gauge for 
various R. 
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FIG. 14: The non-Abelian E electric field profile in the 
Landau gauge obtained with the help of the Wilson loop 
W{R X T — Q X 4:). The solid line denotes the best expo- 
nential fit by the function 1491 1 with the best fit parameters 
given in Table |T] 



Now let us discuss the measurements of the coherence 
lengths. As it was explained above, at least in the LA 
gauge, the operator (or its square-root) is physically 
relevant and may have information about properties of a 
dual Higgs field characterizing the QCD vacuum. Hence 
we expect that the coherence length can be measured 
with the help of ^^(s). In Fig. [THI we show a typical 
example of the Aj^ profile around the string where we 
have adopted the lattice definition (|52|l for the operator 
A^(s). This is very exciting, since the behavior of the 
profile is just what we expect from a profile of a Higgs 
field. 

We plot the R dependence of the coherence lengths and 
the penetration lengths in Fig. 1171 It is very interesting 
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FIG. 17: The coherence lengths of the dimension 2 gluon 
operator and the penetration length in the Landau gauge for 
various R. 
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FIG. 16: The profiles for the dimension 2 gluon operator in 
the Landau gauge around R x T — 6 x 8 Wilson loop. The 
solid line denotes the best exponential fit by the function H49|l 
with the best fit parameters given in Table Q] 



that the values of the coherence lengths are almost the 
same as those for the penetration lengths of non-Abelian 
E electric field. Hence if the observed behavior of the 
correlations is not changed with the increase of R, then 
the type of the vacuum should be fixed to be near the 
border to the type 1 and type 2 as in MA gauge. 



D. Comparison between MA gauge and LA gauge 

In order to consider the gauge-(in)dependence of the 
dual superconductor picture, we show in Fig. ll8l the pen- 
etration length determined in the MA -I- UlLA gauge 
and in the LA gauge. We also plot the coherence length 
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FIG. 18: The penetration lengths of the non-Abelian electric 
field in the Landau gauge and in the MA -I- UlLA gauge for 
various R. 



in Fig. El From these figures, we observe that the coher- 
ence and correlation lengths calculated in different gauges 
coincide with each other. 



V. CONCLUSIONS 

We have observed that 

1. The coherence lengths of the vacuum of the SU(2) 
gluodynamics measured in the Maximal Abelian 
gauge and in the Landau gauge are all consistent 
with each other. 

2. Since the penetration lengths obtained in the 
MA gauge are in agreement with the penetration 
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FIG. 19: The coherence lengths of the dimension 2 gluon 
operator in the Landau gauge and in the MA + UlLA gauge 
for various R. 



lengths calculated in the LA gauge, we conclude 
that the type of the vacuum in both gauges is near 
the border between type 1 and type 2. 

3. The monopole contributions to A'^{s) alone repro- 



duces the full coherence length, although the ab- 
solute value of the correlations is smaller (this last 
fact is quite natural). Therefore the type of the 
vacuum can be determined only from the monopole 
contributions. The observed phenomenon is yet an- 
other example of the monopole dominance. It pro- 
vides an additional support to our expectation that 
the Abelian monopoles are responsible for all non- 
perturbative phenomena related to the confinement 
of color in QCD. 
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